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Exact diagonalization is used to study the quantum states 
of vertically coupled quantum dots in strong magnetic fields. 
We find a new sequence of angular momentum magic numbers 
which are a consequence of the electron correlation in the dou- 
ble dot. The new sequence occurs at low angular momenta 
and changes into the single dot sequence at a critical angu- 
lar momentum determined by the strength of the inter-dot 
electron tunneling. We also propose that the magic numbers 
can be investigated experimentally in vertically coupled dots. 
Because of the generalized Kohn theorem, the far-infrared 
optical absorption spectrum of a single dot is unaffected by 
correlation but the theorem does not hold for two vertically 
coupled dots which have different confining potentials. We 
show that the absorption energy of the double dot should 
exhibit discontinuities at the magnetic fields where the to- 
tal angular momentum changes from one magic number to 
another. 

PACS numbers: 73.2Q.Dx, 71.45. Gm 



Low-dimensional confined electronic systems, such as 
quantum wells, wires and dots, have recently attracted 
much interest because they exhibit dramatic quantum 
effects when they are placed in a strong magnetic field. 
The fractional quantum Hall (FQH) effect which oc- 
curs in a two dimensional sheet of electrons, is probably 
the most spectacular example although there are many 
others. Recently, attention has focused on double-layer 
FQH systems l^]-^ , where the additional degree of free- 
dom (a pseudo-spin that labels the layers) enriches the 
physics. A central issue in these systems is the inter- 
play of electron correlation and inter-layer electron tun- 
neling — because of the competition between these ef- 
fects the quantum Hall state evolves continuously from 
a correlation-dominated (two-component) state down to 
a tunneling-dominated (single-component) state within 
the quantum Hall regime. 

Another recent development is the study of laterally 
confined systems, and quantum dots in particular. A dot 
in the FQH regime contains a few electrons confined on 
a length scale of the order of the magnetic length. As a 
function of the total angular momentum, L, the ground- 
state energy of this system exhibits downward cusps at 
specific L values which are known as the "magic num- 
bers" 0,1). The magic numbers are caused by electron 
correlation, and there are general arguments that relate 
the magic numbers for given numbers of electrons, N^, to 



the symmetry of the wavefunction and the requirement 
that it satisfies Pauli's principle |ll|. For instance, the 
magic numbers, L — 3, 6,9, • • ■, for three spin-polarized 
electrons correspond to triangular spatial correlation that 
minimizes the Coulomb repulsion. 

It is then intriguing to ask what happens if we laterally 
confine a double-layer FQH system to form vertically- 
coupled quantum dots. Because fascinating correlation 
effects are known to occur in double 2D systems (bi- 
layers) and double ID systems (double quantum wires 
|l2| ), we can expect to find interesting phenomena in 
double OD systems (double dots) which are the subject 
of the present work. Specific questions we address are, 
firstly, what will happen to the magic numbers as we 
vary the strength of the tunneling, and secondly, whether 
the magic numbers may become observable in double 
dots. Technically, we believe that the structures con- 
sidered here are within the scope of current fabrication 
technology [ p^ . 

The system wc study is a double dot containing a total 
of three spin-polarized electrons. In both dots the elec- 
tron motion is perfectly two dimensional and the lateral 
confining potential within each layer is assumed to be 
parabolic. The dots are separated in the vertical direc- 
tion with their centers aligned on a common axis. The 
electrons experience both intra- and inter-layer Coulomb 
repulsions in the presence of the inter-layer tunneling. 
We find a new series of magic numbers, which corre- 
spond to ground states dominated by the inter-layer elec- 
tron correlation. Very recently, similar effects have been 
found in l/r^ interaction coupled multiple dots in the 
absence of tunneling [Q. In addition, we propose that 
the magic numbers can manifest themselves in the far- 
infrared optical absorption spectrum. In a single dot this 
is impossible because of the generalized Kohn theorem 
P, P^Jl6| but the theorem does not hold for a double dot 
with different confining potentials. Consequently, we find 
that the absorption energy should exhibit discontinuities 
at the magnetic fields where the total angular momentum 
changes from one magic number to another. 

A vertically-coupled double dot is characterized by the 
strength of the parabolic confinement potential of the 
upper-(lower-)layer, huj^{huj-), the layer separation, d, 
and the strength of the interlayer tunneling (measured 
by AsAs, the energy gap between the symmetric and 
antisymmetric states in the non-interacting system). The 
Hamiltonian, 
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H — Hs + 'Ht + Ti-c, 



comprises the single-electron part, Tig, the tunneling 
term. Tit, and the Coulomb interaction. Tic- We as- 
sume that the magnetic field B is so strong that Landau 
level mixing is negligible and we write the Hamiltonian 
in second-quantized form with a Fock-Darwin basis [ pT[ . 
This gives 
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Here, the index a is used to distinguish the two dots, 
a = , c\^{cea) are creation (annihilation) opera- 
tors and e is the dielectric constant of the host mate- 
rial. The sum over ai, 04 in Tic guarantees that both 
intra- and inter-layer interactions are included. The en- 
ergy of the zeroth Landau-level Fock-Darwin state with 
angular momentum i (> 0) in the ath dot is ega = 
(1 + £)h{uj^/A + ujl)i - ehLUc/2, where iVc = eB/m*c 
is the cyclotron frequency. 

We now estimate typical values of the parameters. To 
obtain approximate values of the confinement energy we 
use a simple electrostatic model in which there is a disk 
of positive charge above and below the two dots, with the 
entire structure sandwiched between two metallic disks. 
This is meant to mimic the electrostatic confinement 
scheme in which a single quantum dot is made by apply- 
ing a modulated gate electrode to a modulation doped 
hetero junction or quantum well For a single dot, 

we have found that the model is able to reproduce the 
confinement energy from an exact solution of the Poisson 
equation to about 20%. To estimate the confinement en- 
ergy of the double dot we take typical device dimensions 
and dopant densities from the work of Boebinger et al 
pO| and Kumar et al jl^ to find that huja is about 2 - 
4 meV. The asymmetry in htOa depends on the offset of 
the two dots from the symmetric configuration and on 
the positions of the disks. It is typically 5-10% for disk 
separations of a few hundred nanometers and dot offsets 
of a few tens of nanometers. Larger asymmetry could 
be achieved by making the structure grossly asymmetric. 
In our calculations we take — 2.0 meV, hcu- = 2.2 
meV. The electrostatic model predicts that the confin- 
ing potential at the center of each dot will be in general 
different. We assume that this could be compensated 
by applying a potential to the entire device. For the 



(1) dot separation and the symmetric-antisymmetric split- 
ting, we take typical values from double layer studies of 
Boebinger et al |^ and Eisenstein et al leading to 
d ~ 20 nm and Asas m the range 0.2 - 0.5 meV. 

The ground-state energy (Fig. |l]) is calculated as a 
function of the total angular momentum, L, by diag- 
onalizing the Hamiltonian in a Slater determinant ba- 
sis at B = lOT for Asas = 0.2 meV (sohd line) and 
Asas — 0.5 meV (broken line). Qualitatively, the be- 
havior shown in the figure is typical for a large range of 
B values although the value of L at which the minimum 
energy occurs depends strongly on B. The magic num- 
bers can be identified from the positions of downward 
cusps. For Asas = 0.2 meV we have a new period of 
two up to L = 3,5,7,9, followed by a period of three, 
L = 9, 12, • ■ •, while for a larger Asas — 0.5 meV the 
period is three throughout, L = 3, 6, 9, 12, • • • as in the 
1^3(^3, iiUi) (jase of a single dot containing three electrons. 

To identify the mechanism for the change of period in 
(4) the magic number for smaller AsaSj we show the charge 
density (the inset to Fig. |l|) and the pair correlation func- 
tion (Fig. H) before (L = 5) and after (L = 12) the change 
in the period sets in. For L = 5 (inset (a)) the density 
against the lateral distance from the center has a peak at 
the center in the lower layer while the density is double- 
peaked in the upper layer. For L = 12 (inset (b)) the 
densities in both layers are double-peaked. We investi- 
gate further by looking at the pair correlation function 
P(r, ro) (Fig. which is defined as the conditional prob- 
ability of finding an electron at position r given that there 
is one at position Tq. The fixed electron is at rg = 16.9 
nm (for L = 5, a) or 23.5 nm (for L = 12, c) in the up- 
per layer where the charge density is a maximum. From 
the figure we can immediately see that the ground state 
electron configuration changes from one dominated by in- 
terlayer correlation to one dominated by intralayer corre- 
lation. For L = 5 the form of the correlation corresponds 
to a triangular " electron molecule" developing across the 
two layers, with one electron at the center of the lower 
layer while the other two are in the upper layer. In con- 
trast, the triangular form develops within each layer for 
i = 12. Similar physics should occur when the two dots 
have the same confining potential (uj^ = w_). 

The change in correlation can be understood by con- 
sidering the energy. As L is decreased the lateral spatial 
extent of the wave function becomes comparable with the 
vertical separation of the layers. When the total angular 
momentum is small enough, the intralayer Coulomb in- 
teraction dominates the interlayer Coulomb interaction, 
so electrons try to avoid each other by developing an 
inter-layer correlation. Although this has to involve mix- 
ing of states in the two dots and costs an energy AsaSi 
the electron correlation still dominates as long as Asas 
is small enough. We believe this is why the new magic 
numbers L = 5,7 appear for smaller Asas- The global 
minimum energy and the angular momentum of the ab- 
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solute ground state depends on the magnetic field. By 
scanning a range of magnetic fields we have found that 
the new magic number states at L — 5 and L = 7 become 
the absolute ground state when B 4T and B ~ 6T re- 
spectively. 

A comparison of our results with the phase diagram 
[Q for the bulk double-layer FQH system is not straight- 
forward. The latter phase diagram is drawn against 
two dimensionless quantities, d/is and AsAs/(e^/e^_B), 
where £b — (ch/eBy^^ is the magnetic length. Be- 
cause of the confining potential the relevant length scale 
for dots becomes the effective magnetic length A with 
= h/m*{Lo1 + icj^y/^. With the parameters we 
have used A = (0.91 - 0.97)^b for B = (5 - 10)T. 
This yields e'^/eX = 14.5 meV for B lOT, so that 
AsAs/(eVeA) = 0.01 - 0.03 for the double dots con- 
sidered here, while the Landau level filling, v, which is 
usually defined as = Ne{Ne — l)/2i for dots, ranges 
from 1/ = 3/5 for L = 5 to z/ = 1/4 for L = 12. It is 
an interesting problem to see how the intra- to inter-dot 
crossover in double dots may be related to the one- to 
two-component crossover in the double layers. 

Now we move on to the far-infrared (FIR) optical ab- 
sorption spectrum. In a single dot with a parabolic con- 
finement potential, the electron-electron interaction does 
not affect the FIR absorption. This follows from the gen- 
eralized Kohn theorem: long- wavelength electromagnetic 
radiations with electric vector E couples to the dot via 
the perturbation Hamiltonian 

N 

H'-^eE-r„ (5) 

which depends only on the center-of-mass coordinate. In 
a single dot with parabolic confinement the Hamiltonian 
separates into the center-of-mass and relative (interac- 
tion) parts and the latter is irrelevant to optical transi- 
tions. In contrast, the separation does not occur in verti- 
cally coupled dots having different confinement energies 
even if both dots have parabolic confinement. 

This means the Coulomb interaction should affect FIR 
absorption spectra. 

To quantify the effect we have calculated the FIR ab- 
sorption spectrum of vertically coupled dots from the ma- 
trix element of the perturbation Hamiltonian, (W), be- 
tween the ground state and all the excited states. Before 
discussing the results, we comment on the applicability 
of this approach to real systems. One important ques- 
tion is the nature of the electric field E. Several authors 
have questioned the relation between the applied electric 
field and the internal electric field in mesoscopic systems 
pl|-p3t with the general conclusion that depolarization 
effects are important. Therefore we would have to cal- 
culate the internal electric field to obtain the absolute 
value of the absorption coefficient. In addition, precise 
calculation of the absorption spectrum would require us 



to take account of other device properties that affect ab- 
sorption, such as finite thickness of the individual dots 
and deviations from a parabolic potential, about which 
scant information is available. We therefore make the 
reasonable assumption that the internal electric field is 
uniform and discuss only the absorption energy and the 
relative intensities of various transitions. This should be 
sufficient for our purpose of demonstrating that the FIR 
absorption of vertically coupled dots is affected by elec- 
tron correlation. 

The results of our calculations (Fig. ^ for Asas = 0.5 
meV shows that the spectrum indeed exhibits a series 
of jumps. In a single dot the FIR absorption has two 
branches: the upper branch for inter-Landau level tran- 
sitions and the lower one for intra-Landau level transi- 
tions. Because we consider only the lowest Landau level 
here we have only calculated the lower branch but we an- 
ticipate that the upper branch will exhibit similar jumps. 
For comparison, the energy of the lower branch for non- 
interacting electrons given by, 

?iWsinglo ^^i^c+ 4^^o)^^^ ~ ^^Wc, (6) 

is also shown in the figure for hcu^ = 2.0 meV (solid line) 
and hiU- = 2.2 meV (broken line). It is clear that the 
coupled dot absorption spectrum is not simple like that 
of a single dot and is split into pieces. This means that 
in the weak magnetic field region, or in a small total an- 
gular momentum region, the center-of-mass and relative 
motions are strongly mixed. 

In particular, the jumps in the absorption energy occur 
at the magnetic fields at which the total angular momen- 
tum changes from one magic number to another. Thus 
the ground state transitions should be directly observ- 
able in the FIR absorption spectrum. The figure also 
shows that the absorption intensity (oc square of the ma- 
trix element) is not monotonic. For Asas = 0-2 meV the 
FIR absorption spectrum is similar, although the jumps 
in energy are smaller than for Asas = 0.5 meV. 

In conclusion, we have found new magic numbers in 
vertically coupled quantum dots and shown that they 
could be probed experimentally. 

We thank Dr N Bruce for providing results on the elec- 
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FIG. 3. FIR absorption spectrum (upper panel) and total 
angular momentum (lower panel) of vertically coupled dots for 
hLj+ = 2.0 meV, hu}- = 2.2 meV, with the layer separation 
d = 20 nm and Asas = 0.5 meV. The position of each filled 
circle gives the energy of the transition while the size of the 
circle represents the absorption intensity. The solid (broken) 
line corresponds to the single-electron absorption spectrum 
for noj± = 2.0 (2.2) meV. 



FIG. 1. Ground-state energy against the total angular mo- 
mentum, L, in vertically-coupled dots with three electrons 
for the strength of the inter-layer tunneling Asas = 0.2 meV 
(solid line) or Asas = 0.5 mcV (broken line). Strength of the 
confinement potential is hui — 2.0 (2.2) meV for the upper 
(lower) layer, and the layer separation is d = 20 nm. Arrows 
indicate the positions of the cusps. The inset shows a cross 
section of the charge density in the upper (solid lines) and 
lower layer (broken lines) against the lateral distance from 
the center of each dot for L = 5 (a) and L = 12 (b). 



FIG. 2. Intralayer (upper panels) and interlayer (lower 
panels) pair correlation functions, P(r, ro), for L = 5, (a,b) or 
L = 12 (c,d). One electron (solid circle) is fixed in the upper 
layer at ro = 16.9 nm {L = 5) or 23.5 nm (L = 12) where 
the charge density has a maximum. The symbol + denotes 
the projection of the solid circle onto the lower layer. An area 
with the linear dimension of 128nm is displayed. The confine- 
ment energies, the layer separation are the same as in Fig.l 
with Asas = 0.2 meV here. 
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